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ABSTRACT 

Analytic expressions of lepton-flavour- and lepton-number-violating decays of charged lep- 
tons are derived in the context of general SU(2)l £g> U(l)y seesaw scenarios that are moti- 
vated by grand unified theories (GUT's) or superstring models, in which left-handed and/or 
right-handed neutral singlets are present. Possible constraints imposed by cosmology and 
low-energy data are briefly discussed. The violation of the decoupling theorem in flavour- 
dependent graphs due to the presence of heavy neutral leptons of Dirac or Majorana nature 
is emphasized. Numerical estimates reveal that the decays r~ — > e _ e _ e + or r~ — > e~ /i~ fi + 
can be as large as ~ 10~ 6 , which may be observed in LEP experiments or other r factories. 
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1 Introduction 



The quest for an understanding of the problem of smallness in mass or masslessness of 
the light known neutrinos, u e , u^, and z/ r , has relied on interesting solutions in the context 
of extended gauge structures of the minimal Standard Model (SM), such as grand unified 
theories, e.g. SO(10) models [[|], or superstring models with an Eq symmetry [0. Among 
the various solutions, the most attractive one, known as the seesaw mechanism, has been 
conceived by the authors in [[J within the framework of SO (10) or left-right symmetric 
models. In these theories, right-handed neutrinos are introduced with the simultaneous 
inclusion of Majorana masses that violate the lepton-number (L) by AL = 2 operators 
in the Yukawa sector. The neutrino-mass spectrum of a simple seesaw model with one 
generation of quarks and leptons consists of two massive Majorana neutrinos, v and N, 
having masses m v ~ m 2 D /mM and mjy — ttim- If the Dirac mass term mp is of the 
order of a typical charged-lepton or quark mass, as dictated by GUT relations [f| , and the 
Majorana-mass scale rriM is sufficiently large, one can then obtain a very light neutrino 
v. The general situation of an interfamily seesaw-type model with a number of weak 
isodoublets and a number n^ of right-handed neutrinos is more involved []5J and will be 
discussed in Section 2. 

If nature keeps to the pathway of a seesaw-type solution, then heavy Majorana neu- 
trinos at the mass scale of TeV may manifest themselves in L- violating processes at high- 
energy ee PJ^] , ep || , and pp colliders |||n| , in possible lepton-flavour- violating decays of 
the Z [ffl| and Higgs particles (H) |12|] or through universality-breaking effects in leptonic 
diagonal Z-boson decays JTH|. Their existence may also influence flIl||T5|| the size of elec- 



troweak oblique parameters [16, IT], tri-gauge boson WW Z- and ZZZ-couplings ||18|, or 



specific Higgs observables considered recently |19|j20|1 . Finally, there are many other places 
scanned by exhaustive combined analyses of charged-current-universality effects in lep- 
tonic 7r decays, neutral-current interactions in neutrino-nucleon scatterings, r-polarization 
asymmetries, neutrino-counting experiments at the CERN Large Electron Positron Collider 
(LEP), etc. [pl] , p2| , in which Majorana neutrinos could also manifest their presence. 

Another possible solution of the neutrino-mass problem has been contemplated in the 



framework of heterotic superstring models || or certain scenarios of SO (10) models |23 
The low-energy limit of such theories extend the SM field content by adding new left-handed 
and right-handed neutral isosinglets, and assuming the absence of AL = 2 operators in the 
Yukawa sector. In a simple one-generation scenario, one obtains three Weyl fermions from 
which one of them is completely massless to all orders of perturbation theory [p^] and the 
other two are degenerate in mass and thus form a heavy Dirac neutrino which has a mass 
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of the order of the isosinglet Dirac mass M. The Dirac mass M simply connects the right- 
handed and left-handed chiral singlets in the Yukawa sector. This solution is particularly 
preferable if the light known neutrinos turn out to be strictly massless. The model could 
straightforwardly be extended to generations without qualitatively changing its features 
regarding neutrino masses. In an n^-generation model, one generally obtains n^ massless 



neutrinos and n^ heavy Dirac neutral fermions p4,25|. This minimal model is invariant 



under the gauge group SU(2)l <8> U(l)y and possesses many attractive features that might 
be summarized in [p6[| . For example, even if the total lepton number is conserved, the 
model does generally violate the separate leptonic quantum numbers and can hence account 
for possible L- and/or CP- violating signals at the Z peak plfl or in other high-energy 
processes |27 . 



In this paper we carefully study the three-body decays of a charged lepton, I, into 
other three charged leptons, which we denote hereafter as V , li, and £2- After detailed 
calculations, we find that the decay amplitude of I — > VIJ2 depends quadratically on the 
mass of the heavy Dirac or Majorana neutrino, which violates explicitly the decoupling 
theorem . Such a nondecoupling behaviour has recently been observed to take place in 
three- generation seesaw-type models, where the effective couplings HIV JL2| and ZW 
show a strong quadratic dependence of the heavy neutrino mass. In past, similar flavour- 
dependent nondecoupling effects have been found in the one-loop amplitude of the decays 
Z —*bs |29| and Z —>■ bb |]30|j3l|] , where the top quark plays the role of heavy neutrinos. 

Among the various decay processes, we find numerically that the decays, r~ — > 
e~e~e + and r~ — > e~fi~fi + [or complementary the decays r~ — > fi~fi~fi + and t~ — > 
fi~e~e + ], have the biggest opportunity to be detected at the present or future LEP data. 
Furthermore, we analyze the effect of genuine Majorana- neutrino contributions to these 
decays. 

The present work is organized as follows: In Section 2, we give a brief description of 
the basic low-energy structure of the seesaw-type models mentioned above. In Section 3, we 
discuss general constraints that should be imposed on these models. Analytically, Section 
3.1 considers possible constraints based on the assumption that the model should generate 
a sufficiently large lepton asymmetry via the out-of-equilibrium L-violating decays of a 
heavy Majorana neutrino which can be converted later on into the observed baryon-number 
(B) asymmetry in the universe (BAU) due to the sphaleron interactions. In section 3.2, 
stringent constraints of possible non-SM mixings are derived by a global analysis of all 
existing low-energy data. Also, bounds that may be obtained by the non-observation of 
leptonic non-diagonal Z-boson decays at LEP are discussed in Section 3.3. In Section 4, we 
analytically calculate the branching ratios of the photonic decays of a lepton (I), I — > I'j, 
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and the three-body decay modes of the type / — > I'lifa in the context of the models discussed 
in Section 2. Numerical predictions and discussion of these lepton-flavour- violating decays 
are summarized in Section 5. We draw our conclusions in Section 6. 



2 Theoretical models 

In this section, we will give a short description of the basic low-energy structure of the 
two most popular extensions of the SM that can naturally account for very light or strictly 
massless neutrinos. The field content of these models, which could also be motivated by 
heterotic superstring models |2| or certain S'O(10)-GUTs |||23|], is free of anomalies |26| . 



These two scenarios are: (i) the interfamily seesaw-type model realized in the SM with 
right-handed neutrinos HHH and (ii) the SM with left-handed and right-handed neutral 
singlets |0f2§[25|. 



(i) The SM with right-handed neutrinos. In general, the interfamily seesaw- type 
model, being invariant under the SM gauge group, represents one of the most natural 
framework to predict heavy Majorana neutrinos. Such a model is obtained by introducing 
a number n# of right-handed neutrinos, z/^, in the SM [in addition to n^ left-handed ones 
and allowing simultaneously the presence of AL = 2 operators. The Yukawa sector 
containing the neutrino masses is then written down as 

1 / j.c 

-JCXr = -fi/9, vl c )M»\ \ | + B.C., (2.1) 




where the (n^ + iir) x (uq + n^)-dimensional neutrino-mass matrix 

M" = ( ° T mD ) . (2.2) 
y mjj m M J 

The matrix M v can always be diagonalized by a unitary matrix U u of the same dimen- 
sionality with the neutrino-mass matrix (i.e., U vT M V U V = M v ). One then gets n^ + njj 
physical Majorana neutrinos rij through the unitary transformations 

v \ iic+nji / OC \ n G +n fl 

Jc = E U£ n Lj and ^ = ^ U» n Ry (2.3) 

The first n^ neutral states, z/j (= ni for i — 1, . . . , io.q), are identified with the known 
light neutrinos (i.e., n^ = 3), while the remaining np mass eigenstates, Nj (= nj +riG for 
j = 1, . . . , tir), represent heavy Majorana neutrinos which are novel particles predicted by 
the model. The quark sector of such an extension can completely be described by the SM. 
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It is important to notice that the general matrix M v of Eq. (2.2) takes the known 
seesaw form || in case 3> mjy. Nevertheless, this hierarchical scheme can drastically 
be relaxed in a two-family seesaw-type model without contradicting experimental upper 
bounds on light-neutrino masses @J^J§]. The light-heavy neutrino mixings of such scenar- 
ios, s^, can, in principle, be scaled as ~ mo/mjf rather than ~ Jm u JmN as usually 



derived in a one-family seesaw scenario HJ32|. In other words, high Dirac mass terms are 



allowed to be present in M y and only the ratio mo/m^ (~ s 2) & e ^ s limited by a global 
analysis of low-energy and LEP observables. The latter advocates our treatment of origi- 
nally considering the mixings s2 and heavy neutrinos masses m^. as free phenomeno logical 
parameters, being subject later on to the constraints that will be discussed in Section 3. 

Adopting the conventions of Ref. ||, the interactions of the Majorana neutrinos, rii, 
and charged leptons, k, with the gauge bosons, W ± and Z, and the unphysical Goldstone 
bosons, G ± and G° (in the Feynman-'t Hooft gauge), are correspondingly obtained by the 
Lagrangians: 



C 



w 

hit 



n G n G +n fi 



£ r Ih jIr.V,. H.c. 



fju 



i=l j=l 



n G +n H 



4c 



n.-, 



(2.4) 
(2.5) 



and 



* — • , , 



■hit 



9u 



V2M } 

l 9w ^yO 



n G n G +n R 
G ~ zl /J B hj k 



U i=l j=l 

n G +n fl 



rij + if.c. 



mi 



4M, 



(2.6) 



(2.7) 



where g> w is the weak coupling constant, c 2 w = 1 — s 2 w = M^/M^, P L (P R ) = (1 + (— )7s)/2, 
and Wj denotes all the physical neutrino masses. In Eqs. (2.4)-(2.7), B and C are n G x 
(n R + n G )- and (n G + n R ) x (ji R + n G )-dimensional matrices, respectively, which are defined 

as 

n G n G 
Bid =12KkU% and Qj = £ UkPiS> ( 2 - 8 ) 

fc=l h=l 

where V 1 is the leptonic Cabbibo-Kobayashi-Maskawa (CKM) matrix. 

Note that the flavour-mixing matrices B and C satisfy a number of identities, which 
are derived just by using the information of SU(2)l <8> U(l)y invariance of C Y - These 
identities, which are forced by the renormalizability of the interfamily seesaw-type model, 
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can be summarized as |12|,p| 



2J B hkB* hk 



k=l 



k=l 

n G +n fl n G +n fl 



n G +n fl n G 

k=l k=l 
n G +n fl , 



^ m k C ik C jk = 0, ^ m k B lk C* ki = 0, ^ m k B hk B, 



l 2 k 



0. 



a,-, (2.9) 

(2.10) 



fc=i 



fc=i 



fc=i 



Consequently, our theoretical analysis should be regarded to be independent of the weak- 



basis structure of possible neutrino- mass-matrix ansatze | 33fl . It is now instructive to 
re-express the Z-boson coupling to the Majorana neutrinos, rii, as follows: 



'-'mt 



4c, 



i,3=l 



(2-11) 



One can thus remark that the coupling Zn^nj is generally flavour non-diagonal and has 
both chiralities in this minimal model. 

(ii) The SM with left-handed and right-handed neutral singlets. Another attractive 
scenario predicting for the light neutrinos to be strictly massless serves an extension of the 
SM, in which left-handed neutral singlets, Su, in addition to the right-handed neutrinos, 
2/£ , are introduced. Furthermore, we assume that AL = 2 interactions are absent from 
the model, and the number of right-handed neutrinos, n#, equals the number of the singlet 
fields Sn- After the spontaneous break-down of the SM gauge symmetry, the Yukawa 
sector relevant for the neutrino masses reads 



-a 



Y 



y~Li v r > 



S L )M U 



( 



V 



oc 



\ 



'R 
C 



+ H.c, 



(2-12) 



where the {uq + 2n#) x (% + 2n#) neutrino-mass matrix is given by 



( o 



m 



V o 



m D 


M 





M T 




\ 



(2-13) 



It is easy to see that the neutrino matrix M. u conserves the total lepton number L by 
assigning the lepton numbers to the neutrino fields: L^v®) = L(v R ) = L(Sl) = 1. Since 
the rank of the neutrino mass matrix in Eq. (2.13) is 2n#, Ai v has n^ zero eigenvalues. 
These n G massless eigenstates should clearly describe the ordinary light neutrinos, v e , 
and v T p|,[23|] . The remaining 2nR Weyl fermions are degenerate in pairs due to the fact that 
L is conserved and so form njj heavy Dirac neutrinos. In general, this viable seesaw-type 



model can have large Dirac components in J\A U and only the ratio mn/M 



~ s. 



will again 
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be constrained for M~ 100 GeV (see also our discussion in Section 3). A nice feature of the 



model is that the individual leptonic quantum numbers may be violated |2^|25|J34| , even if 
L is conserved. The charged-current Lagrangian can be obtained by Eq. (2.4), while the 
neutral-current interaction is given by 

n G +n R 

Cf nt = - f^Z» £ o.Cr^U'j. (2.14) 

The matrices B and C for this specific model obey the sum rules in Eq. (2.9), but not the 
identities of Eq. (2.10). 

At this stage, we must comment on the difference between the Lagrangians (2.5) 
and (2.14). Since Eq. (2.5) describe Majorana neutrinos contrary to the Lagrangian (2.14) 
where the massive neutrinos are Dirac, the strength of the Zfiifii coupling for identical 
Majorana fermions is two times larger than the one which may naively be red off from Cf nt 
in Eqs. (2.5) and (2.11). The off-diagonal coupling ZfiiTij (with n« ^ rij) is again enhanced 
by a factor of two, since the charge-conjugate interaction Znjrii will equally contribute to 
the coupling of the Z-boson to Majorana neutrinos. In our forthcoming calculations, we 
have taken into account all these theoretical differences in treating Majorana and Dirac 
fields. In fact, we find that taking formally the limit C*j — > but keeping CV, 7^ in 
Eq. (2.11) and considering the afore-mentioned statistical Majorana factors is sufficient to 
recover the model with additional left-handed neutral singlets. 

To make life easier, we ultimately make the following reasonable assumptions: To a 
good approximation, we assume that possible novel particles related to the above unified 
theories, such as leptoquarks or extra charged and neutral gauge bosons, W% and 
Z' [p7| , ^2| , are sufficiently heavy so as to decouple completely from the low-energy processes 
discussed in Sections 3.3 and 4. For obvious reasons, possible singlet and triplet Majoron 



fields [p8| , |39| , ^0| , ^ll are considered to couple very weakly to matter so that we can safely 
ignore them in our considerations. 



3 General constraints on the models 

3.1 Cosmological constraints 

Unified theories based on the gauge group 50(10) or Eq can naturally accommodate 
right-handed neutrinos in addition to quarks and leptons of the SM. In such theories, the 
Majorana mass tum can directly be related to the B — L scale of a local symmetry which 
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is assumed to be spontaneously broken. It is therefore evident that the mass of heavy 
Majorana neutrinos will be determined by the scale of B — L breaking. Moreover, out- 
of-equilibrium lepton-number-violating decays of heavy Majorana neutrinos can generate 
a non-zero L in the universe through the L- violating interactions of Eqs. (2.4)-(2.7). 
This excess in L can be converted into a B asymmetry of the universe via the (B + 
L)-violating sphaleron interactions, which are in thermal equilibrium above the critical 
temperature of the electroweak phase transition [fi"3| , fi4}| . Many studies have recently been 
devoted to constrain the (B — L)-violating mass scale by making use of the drastic out-of- 
equilibrium condition for the AL = 2 operators, and so to derive a lower mass bound on 
the heavy Majorana neutrinos [^5|j4^j47|j4^j4^|] . For example, in conceivable scenarios 
predicting heavy Majorana neutrinos with = 1 — 10 TeV could naturally account for 
the observed BAU. Subsequently, it was argued [ffl that the rriN lower bound of ~ 1 TeV 
was considerably underestimated and a lower bound on%> 10 5 TeV should be imposed 
in a two-generation scenario of right-handed neutrinos with large interfamily mixings so 
as to be compatible with the existing BAU. This would obviously imply that probing 
Majorana-neutrino physics at collider energies may not be phenomenologically interesting. 

The latter observation can indeed be valid in a two-generation-mixing model with 
two right-handed neutrinos. In general, in three-generation models with lepton-flavour 
mixings, a careful inspection of chemical potentials has shown that the stringent mass 
bound of heavy Majorana neutrinos mentioned above can be weakened dramatically and 
is quite model dependent f48],f49|. In particular, it is sufficient that in equilibrium state 
one individual lepton number, e.g. L^, is conserved in order to generate the BAU via the 
sphaleron interactions, even if nonvanishing operators with AL/ t ^ AL^ were in thermal 
equilibrium |j48| . The reason is that sphalerons generally conserve the quantum numbers 
B/3 — L[ t [050 an d thus preserve any BAU generated by an excess, e.g., in L M , from 



being washed out. Similar conclusions have been drawn in Ref. p*f|. 

A viable scenario of heavy Majorana neutrinos with masses ~ TeV can easily be 
realized in the SM with nfj = 4. If the BAU is to be generated through an excess in 
muonic number density, this asymmetry in L M can be achieved by considering a neutrino 
mass matrix, M u , similar to the CP- violating scenario given in ||19|| . This scenario contains 
one left-handed neutrino, z/£ 1; to which, for the case at hand, a muonic quantum number 
should be assigned, and two right-handed neutrinos, denoted as z/^ 3 and z/^ 4 . The explicit 
form of M v is then given by |19 

/ a b \ 



a A 
V b 




B J 



(3.1) 



S 



where a and b are in general complex numbers, and A and B can be chosen to be real. 
Out-of-equilibrium conditions for generating a sufficiently large asymmetry in L M , which 
can give rise to the established BAU, lead to the stringent lower bounds on the masses of 
the corresponding physical heavy neutrinos as consistently obtained by |I6] . However, 
the remaining e- and r-lepton families can strongly mix each other via two additional 
right-handed neutrinos, e.g., v° R1 and u R2 , and form an individual 4 x 4 seesaw-type matrix. 
Operators AL e ^ and AL T ^ are now allowed to be in thermal equilibrium provided 
that A(L e — L /t ) = and A(L r — L M ) = 0. The latter condition is automatically satisfied 
due to the construction of this specific scenario with iir = 4. As a consequence, the severe 
lower mass bounds on the physical heavy neutrinos N% and N 2 can be evaded completely. 
A similar analysis in a SM with n# = 3 is more involved due to the flavour-mixing effects 
in the neutrino-mass matrix and can be given elsewhere. 



3.2 Low-energy constraints 

There exists a great number of low-energy experiments that could set upper bounds 



on possible non-SM couplings pl| . The most significant experimental tests giving strin- 
gent constraints turn out to be the neutrino counting at the Z peak, the precise measure- 
ment of the muon width fi — > ev e v^, charged-current universality effects on the observable 
r(7r — > ez/)/r(7r — > fiu), non-universality effects on B(t — > evv)jB{r — > [ivv\ and other 
nuclear physics effects and experiments. All these constraints, which are derived by the 
low-energy data mentioned above, depend crucially on the gauge structure of the model 
under discussion. For example, assuming the supersymmetric (SUSY) nature of the E 6 
models or SUSY-S'O(IO) unified theories [50,pl|, and i?-parity invariance, the neutralino 



state could then be the lightest supersymmetric particle (LSP) which is stable. If the mass 
of the LSP is assumed to be in the vicinity of Mz/2, then an additional invisible decay 
channel for the Z boson will open kinematically and neutrino-counting limits imposed on 
the couplings Zi/iUj may not be applicable. Furthermore, an analysis of decays of the type 



Z — > Nu, which have been considered in f27j, suggests that 100 GeV for (s^) 2 ~ 0.01. 



Thus, identifying the non-SM-mixing angles {s2) 2 of Ref. |2T[] as 



= £IW> (3-2) 



and in view of the discussion given above, one may tolerate the following upper limits |21 



s v lf < 0.015, (s v £) 2 < 0.070, and (s^) 2 < 1. 10" 9 . (3.3) 
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In Eq. (3.3), the tight upper bound on s^ 1 represents that the muonic quantum number 
is practically conserved in thermal equilibrium. Note that, without loss of generality, one 
could equally interchange the upper limit on (s v £) 2 with that of (s^) 2 . To be precise, we 
will assume (s^ 1 ) 2 ~ 0, in what follows. 

Another limitation to the parameters of our model comes from the requirement of the 
validity of perturbative unitarity that can be violated in the limit of large heavy-neutrino 
masses. A qualitative estimate for the latter may be obtained by requiring that the total 
widths, T Ni , and masses of neutrino fields iVj satisfy the inequality 

^ < I (3.4) 

m Ni 2 

The total widths of the heavy neutrinos, T^, can be written down as a sum over all possible 
decay channels ||, i.e. 

r Ni = E r (^ - l t WT ) + £ (iW - VjZ) + iW - Vj H)). (3.5) 

In the limit of ^> Mw, Mz, Mjj, the above expression simplifies to 

r ^ = a J2 m N, \^NiNi\ 2 , (3-6) 



4M 2 , 



with a w = <7^/47r. 



3.3 Constraints from leptonic Z-boson decays 

Aside from low-energy constraints discussed in Section 3.2, many extensions of the 
SM derived by unified theories may give rise to lepton-flavour-violating decays of the Z 
boson [ 52 J53]|LT] . In particular, it has been found in |Tl[ that the non-observation of such 



non-SM signals at LEP may impose combined bounds both on heavy neutrino masses 
and mixings (s^) 2 . The reason is that the amplitude of such a decay depends quadratically 
on the heavy neutrino mass, leading to measurable rates. In a self-explanatory way, the 
amplitude of the decay Z — > IV may generally be parametrized as 

T(Z-H') = |^ Ti 4 %7,(1 - T5K (3.7) 

where a w = g^/^n and the form factor JF|', which is induced by the Feynman graphs of 
Fig. 2 at the one-loop electroweak order, is given in Appendix A. The branching ratio of 
this decay mode is obtained by 
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where V z = 2.49 GeV is the experimental value of the total width of the Z boson |56[. We 
postpone the numerical discussion of possible constraints that might arise due to lepton- 
flavour-violating decays of the Z boson in Section 5. 



4 Flavour-violating decays of charged leptons 

In Sections 4.1 and 4.2, we will theoretically analyze the possibility of lepton flavour 
nonconservation in decays of the form I — > £'7 and I — > respectively. As mentioned in 
the Introduction, /, l x and l 2 denote usual charged leptons, i.e. the e, \i and r leptons. 



4.1 The decay I — > I'j 



In the framework of the minimal class of models discussed in Section 2, heavy Majo- 
rana or Dirac neutrinos can give rise to the decay I —>■ l /r y. The Feynman graphs responsible 
for such a decay are shown in Fig. 1. Applying electromagnetic gauge invariance to the 
decay amplitude l{p) — > l'{p')^{q), where the photon, 7, may be off-mass shell, yields [|55l ,|57 



GOi w _ 



1 w 



n G +n R 

£ 5f^F 7 (A,)(g 2 7 M 
i=i 



-75) 



- B^BuiG^Xi) ia^ v q u [m v (l - 75) + mi(l + 



75 



i=l 



Ul, 



(4.1) 



where Aj = mf/M^, q = p — p' denotes the outgoing momentum of the photon, and the 
form factors F 1 and Gy are given in Appendix B. It is now straightforward to calculate the 
branching ratio of I — > I'j 



b(i - r 7 ) 



ajsl m\ mi w |2 
2567r 2 M^r, 1 t 



(4.2) 



where Ti is the total width of the decaying lepton I, while Gj: in Eq. (4.2) represents a 
composite form factor defined in Appendix B. Specifically, for the total width of the r 
lepton, we use the experimental value, T T = 2.1581 1CT 12 GeV [56 1, whereas the muon 
total width may be obtained by |53| 

,2\ r ,25 



1927T 3 



1 - 



m. 



1 



2tt V 4 



7T 



(4.3) 



where a em = e 2 /47r. The muon total width given in Eq. (4.3) is in excellent agreement 
with the experimental value reported in [5B|. 
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The experimental upper bounds arising from the non-observation of decays of the 
type I — > Z'7 are 

B(r -> e 7 ) < 1.2 1(T 4 , B(r -> ^7) < 4.2 1(T 6 , 5(/i -> e7) < 4.9 KT 11 , (4.4) 

at 90% CL. Using the values for the mixing angles (s^) 2 of Eq. (3.3), one easily finds that 
the photonic decays involving muons are extremely suppressed in our minimal scenarios. 
Furthermore, the theoretical prediction B(r — ► e 7 )~ 1CT 7 shows that photonic decays of a 
r lepton may not be the most favourable place to probe heavy neutrino physics. 

4.2 Three-body leptonic decays I — > Vl\l2 

In a three-generation model the decaying charged lepton I will either be a muon or a 
t lepton. There are seven possible decays of the generic form I —>■ l'lil 2 

a. r~ — > 

b. r~ — > fi~{i~e + , 



C. T — ► 6 fi 

d. r~ — > t~ t~ yU 4 



e. r — > e /x e , 

f. t~ — > e~e~e + , 

g. /i~ — > e~e~e + . (4.5) 

To facilitate our computational task, we divide the decays in Eq. (4.5) into three categories 
according to the leptonic flavours in the final state: Category (i) contains all the decays 
where V 7^ l 2 and l\ = l 2 or V = l 2 and l\ 7^ l 2 {i.e. the decays (c) and (e)). Category (ii) 
comprises all the decays where V = l\ = l 2 (i.e. the decays (a), (f) and (g)). And lastly, 
all the decays with final leptons having I' 7^ l 2 , l\ 7^ l 2 belong to the category (iii) (i.e. the 
decays (b) and (d)). 

The transition amplitude of the decay l(p) — > I' (p')h(pi)h(p2) receives contributions 
from 7- and Z-mediated graphs shown in Figs. 1 and 2, respectively, and box diagrams 
given in Fig. 3. These three different amplitudes are conveniently written down as follows: 



T 7 (l^l'hl 2 ) = --^S hh £ B*5 w %7X%^7(Ai)(7M-^f)(l-T 5 ) 
^ lvl w i=i Q 

-zG 1 (X i )a flu — (mi(l + 75) + m v (l - 75)) u h (4.6) 
q z J 



2 c 2 n G +n fl 
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x5 hh J2 B u B i'i [^Fz(Xi) + C i3 H z {\ u Xj) + C*p z {\, \)\ ,(4.7) 

„2 



T Bo x{l -> Z'Zife) = - 16i J 2 "Z'T/X 1 - 7s)^ ^7^(1 - 7 5 )v, a 

+Bi'iBhiB*jB* a j Gbox{K, Aj) , (4.8) 

where g = pi + in addition to the photonic form factors _F 7 and G 7 in Eq. (4.6), the 
form factors F^, ii^, G^, F Box , and Gboz are given in Appendix B. Note that the term 
proportional to G 7 in Eq. (4.6) contains a non-local interaction which is singular in the 
limit q 2 — > 0. 

Following the classification mentioned above, the branching ratio for all decays be- 
longing to the first category is found to be [] 

bq- -> r i7i+ , v ^ i 2 , h = i 2 ) = a ™ - ^i- — 

V 1 2 ' T 2 ' 1 2; 24576tt 3 T z 

x { I Ff ^ + if - 2sUF§' - if) | 2 + 44 I if - if | 2 
+ 84 Re[(if + Ff^)G«*\ - 324 He[(if - if )<f *] 

+ 324IG?'! 2 [In 4 -3]}, (4.9) 



where if, Cf, if, and if j^ 2 are composite form factors defined explicitly in Appendix 

B. The branching ratios referring to the categories (ii) and (iii) are correspondingly given 
by 

K 1 2 ' ; 24576tt 3 Ti 

v/o I 1 p'Wl'l _l_ 7?Hl o„2 / pKi J?lh\ 1 2 1 /| „4 I jplli zpllj |2 

+ 164Re[(if + \F l ^ h )G 1 ^) - 484Re[(F^-Ff)G^*] 

+ 3241^ | 2 P^^-tI} ( 4 - 10 ) 



*In our calculations we have used a notation similar to the authors in [ p4j . However, our branching-ratio 
expressions (4.9) and (4.10) are at variance with their results. 



13 



and 

b<t + h ,e ,« - ASS I i 2 ■ <*"> 

Equations (4.9) and (4.10) contain a non-local interaction in terms oc G"' and C" 1 , which 
is discussed in detail in Appendix C. In Eq. (4.10), one has to take into account statistical 
symmetrization factors for the two identical final leptons (i.e. I' = li), as well as additional 
Feynman graphs resulting from the interchange of the lepton I' with h. The set of decays 
in (iii) can only be induced by the box graphs shown in Fig. 3. The amplitude of such 
a class of decays (i.e. decays (b) and (d)) is always proportional to (s 1 ^) 2 s v £ and the 
corresponding branching ratios are hence expected to be vanishingly small even if one uses 
the upper value of s v £ in Eq. (3.3). For reasons of mere academic interest, we simply note 
that B(t~ — > e~e~/z + ), B(t~ — > yU _ /i~e + ) ~ 10 -12 . As a consequence, we find that the 
decays (c) and (f) in Eq. (4.5) deserve the biggest attention and will hence be analyzed 
numerically in the next section. 



5 Numerical evaluation and discussion 

We will now investigate the phenomenological impact of the two seesaw-type models 
outlined in Section 2. In order to pin down numerical predictions, we will, for definiteness, 
assume an extension of the SM by two right-handed neutrinos. The neutrino mass spectrum 
of such a model consists of three light Majorana neutrinos which have been identified 
with the three known neutrinos, v e , v^, and v T , and two heavy ones denoted by N\ and 
N 2 . As already mentioned in Section 2, the seesaw- type extension of the SM with one 
left-handed and one right-handed chiral singlets can effectively be recovered by the SM 
with two right-handed neutrinos when taking the degenerate mass limit for the two heavy 
Majorana neutrinos. It is therefore obvious that branching-ratio results for the SM with 
one left-handed and one right-handed neutral singlets can be red off from the SM with two 
right-handed neutrinos in the specific case m Nl = m^ 2 = m N . 

Apart from the two heavy Majorana neutrino masses which are free parameters of 
the theory, the model contains numerous mixing angles, Bu and Cy, for which the only 
restriction comes from a low-energy analysis as discussed in Sections 3.2 and 3.3. However, 
in our minimal model with two right-handed neutrinos one can derive, with the help of the 
identities in Eqs. (2.9) and (2.10), the useful relations 

p ^~ ^ S £ 

Bin! = I , B m = - , (5.1) 
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where p = m 2 N2 /m 2 Ni . The mixings C^n^ can also be obtained by employing Eq. (2.9). In 
this way, one gets 



pl/2 n G j n G 

i 7 „i/ 2 E( s l) 2 ; Civ 2 w 2 = 1/2 ^K0 2 > 



C^i^a = -Cn 2 Ni = -f- 777 XX s l) • ( 5 ' 2 ) 

l+pV2 ^ 

Evidently, our minimal scenario depends only on the masses of the heavy Majorana neu- 
trinos, rriN 1 and mjv 2 (or equivalently on and p), and the mixing angles (s^) 2 , which 
are directly constrained by a global analysis of low-energy data. 

In our illustrative model, with the help of Eq. (5.2), we can easily obtain the maximal 
heavy neutrino mass allowed by perturbative unitarity. Satisfying Eq. (3.4) for both heavy 
neutrinos N\ and N2, one gets the global relation 

-i 



2M 2 W 1+p- 1 / 2 



n G 



E( 

i=i 



(5.3) 



with p > 1. Condition (5.3) has thoroughly been used in our numerical estimates to impose 
an upper bound on m,N 12 - 

For reasons mentioned in Section 4.2, we present the branching ratios for the leptonic 
decays — > e~e~e + and r~ — > e~p~ p + in Fig. 4. To gauge to which extend our minimal 



model can predict measurable rates, we have first assumed the maximally allowed values [^1 
for (s v £f = 0.07 and (s u L e ) 2 = 0.015 {{s u £f ~ 0) given in Eq. (3.3). From Fig. 4 we find 
the encouraging branching ratios 

B(t~ -> e- e -e + ) ~ 2. 10" 6 and B(t~ -> e»+) ~ 1. 10~ 6 . (5.4) 

The present experimental upper limits on these decays are given by fl5"6"| 

5(r~ -> e _ e~e+) < 1.3 10" 5 and B(r" -> e» + ) < 1.9 10~ 5 , CL = 90%. (5.5) 

Even if we assume smaller values for the mixing angles, (s v £) 2 = 0.035 and (s^ ) 2 = 0.01 
(( s zf) 2 = O)' the lepton-flavour-violating decays of the r lepton can still be significant. 
From Fig. 5 one has that 

B(t~ ->■ e~e~e + ) ~ 5. 10~ 7 and B(r~ -> e"p"p + ) ~ 3. 10~ 7 , (5.6) 

and the possibility of observing such decays at future LEP experiments appears feasible. 
Note that the branching ratio increases with the heavy neutrino mass to the fourth power 
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and hence allows to reach measurable values. To demonstrate this fact, we have just ne- 
glected contributions of seemingly suppressed terms C((s^) 4 ) in the transition elements 
and found a reduction of our numerical values up to ~ 10~ 2 . In the low-mass range of 
heavy neutrinos (i.e. for < 200 GeV) the difference between the two distinct compu- 
tations is quite small and consistent with results obtained in ||34j| . In the high- mass regime, 
however, the situation changes drastically (see also Figs. 4 and 5), since in the transition 
amplitude, terms proportional to (s2) 2 depend logarithmically on the heavy neutrino mass 
mjv, i.e. \n(m 2 N / M^y) ■, while terms of 0((s u l) A ) show a strong quadratic dependence in the 
heavy neutrino mass, i.e. mj^/M^. 

Fig. 6 represents genuine Majorana- neutrino quantum effects, since we have computed 
the branching ratios as a function of the ratio m^/m^ for the selective values of = 
200 GeV and 500 GeV. Although the most stringent constraints on the heavy Majorana 
neutrino masses result from Eq. (5.3), it is, however, important to notice that for lower 
neutrino masses the maximum of B(r~ — > eTe~e + ) and B(r~ — > e~/i~/i + ) is not given 
by the degenerate case where p = 1. In fact, if m^/m^ ~ 3 the branching ratios show 
up a maximum which can be up to two times bigger than the case where both the heavy 
neutrinos, Ni and N 2 , are degenerate. We have thus found that for = 500 GeV and 
m N2 ~ 1.5 TeV, 

B(t~ -> e"e-e + ) ~ 2. 10~ 8 and B(t~ -> e»+) ~ 1.5 HT 8 . (5.7) 

Such effects might be accessible at r factories if one assumes an upgrade in the luminosity 
of the LEP collider by a factor of 10. 

In the following, we will discuss possible constraints that might arise from lepton- 
flavour-violating decays of the Z boson. Since we always assume that (s u £) 2 = 0, we will 
focus our analysis on the decays Z — > e~r + + e + r~. Within the perturbatively allowed 
range of heavy neutrino masses as determined by Eq. (5.3), Fig. 7 gives 



B(Z - 


- + i + — \ 
-> e r + e t ) ~ 


4.0 10" 6 , 


for 




) 2 = 0.070, 




) 2 = 0.015, 


B(Z - 


- + i + — \ < 
-> e t + e r ) ~ 


1.1 10" 6 , 


for 




) 2 = 0.035, 




) 2 = 0.010, 


B(Z - 


-> e _ r + + e + r~) ~ 


6.0 10" 7 , 


for 




f = 0.020, 




f = 0.010 



In our numerical estimates, we have used values for (s^) 2 compatible with the updated 
upper bounds given in Eq. (3.3). Although all the branching ratios in Eq. (5.8) could be 
detected at future LEP data, they cannot impose any severe constraints on the r decays 
into three charged leptons for the present analysis. The experimental sensitivity at LEP is 
currently given by jpBfl 

B(Z -> e~r + + e + T~) < 1.3 10~ 5 , CL = 95%. (5.9) 
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Here, some comments are in order. In Fig. 7, the branching ratios for the three different 
mixing-angle sets in the order stated in Eq. (5.8) show a minimum at the positions = 
700, 900 and 1200 GeV, respectively. The reason is that C((s^) 2 ) and C((s^) 4 ) terms 
of T l z in Eq. (3.7) cancel each other and the whole transition amplitude becomes pure 
absorptive. In the range of very heavy neutrinos, terms proportional to (s^) 4 will dominate 
in the amplitude for the same reasons mentioned above. The effect of such a dynamical 
cancellation of the dispersive part of the amplitude could be shown up as a difference 
between the charge-conjugate decay modes of Z — > e~r + and Z — > e + r~, leading to 
sizeable CP-violating effects . 

In Fig. 8, we display genuine Majorana- neutrino virtual effects by examining the 
behaviour of the branching ratio as a function of the quantity m^/m^ for rather modest 
values of m^ 1 . Here, the situation is more involved and depends strongly on the value 
of rriNi we choose. The fact that the amplitude could be dominated by (s^) 2 terms for 
relatively light heavy Majorana neutrinos {i.e. rriN 1 < 400 GeV) or by (s^) 4 terms for 
larger values of plays a crucial role for the shape of the different lines drawn in Fig. 8. 
The common feature is, however, that the case where both the heavy Majorana neutrinos, 
N\ and N2, have the same mass does not again correspond to the situation yielding the 
biggest branching-ratio value. 

Finally, r leptons can also decay hadronically via the channels: r — > lit], r — > Z.;7r°, 
etc. ||34j| . However, the present experimental sensitivity to these decays seems to be rather 
weak [56|] so as to set constraints on our analysis. For example, B(t — > en ) < 1.4 10~ 4 , at 
CL= 90%. 



6 Conclusions 

We have explicitly shown that seesaw-type extensions of the minimal SM, which 
naturally contain left-handed and/or right-handed weak isosinglets, can favourably account 
for sizeable branching ratios of r decays into three charged leptons that can be as large as 
~ 10 -6 . Using updated constraints for the mixings (s^) 2 , we have found that our numerical 
estimates of B(t — > eee) and B(t — > e^ifi) are in qualitative agreement with those obtained 
in [34 for < 200 GeV. However, our branching-ratio values can be up to 100 times 



larger than the results reported in |34| when the heavy neutrinos have TeV masses. The 



reason is that the flavour-violating decays of the r lepton show a strong quadric mass 
dependence of the heavy neutrino mass in a complete calculation, which gives a unique 
chance for such decays to be seen at LEP or planned collider machines. 
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Apart from general constraints that our minimal models should satisfy and have been 
taken into account, we have found that B(Z — > e~r + + e + r~)~ 4. 1CT 6 within the range 
allowed by perturbative unitarity. The latter do not yet impose any stringent constraints 
on the phenomenological parameters of the theory. Moreover, heavy Majorana neutrinos 
introduce a different behaviour in the transition amplitude via loop effects as compared 
to heavy Dirac ones. For example, Fig. 6 shows that an appreciably large mass difference 
between the two heavy Majorana neutrinos Ni and N 2 (i.e. m^/m^ ~ 3) will give rise 
to an enhancement of a factor of two to the corresponding branching-ratio value obtained 
for = TTi jv 2 • 
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A Loop integrals of leptonic Z-boson decays 



After computing the Feynman graphs shown in Fig. 2, we find that the analytic 
expression of the form factor T l z defined in Eq. (3.7) can be cast into the form [|TT[ 



n G +n R 



•J 7 



i,j=i { 



--(l-2s' w )X i [2L 1 (X i ) + 



I(Xi) - 3c^Li(Aj) - s 2 w \il(\i) 

3 3 (A l + 2)A l lnA i 



2 1 - A,- 



+C y (^L 2 (\, Xj) - ^Xz[K x {\, X 3 ) - K 2 (\, X 3 ) + K(X h A,)] - ^X^K^Xi, X L 



+Cy\Xj l-K^Xi, Xj) + -X z K(Xi, A,) - -L 2 (Xi, A/ 



(A.l) 



where Aj = mf/M^, Xz = M^/M^-, and the definition of the loop integrals, /, I, Li, K%, 
K 2 , K, and L 2 , may be found in [fL3|j . The analytic expressions of these loop integrals are 
listed below 



Lx{Xi) 



1 rl 



dxdyy 



1 



o Jo Bx(Xi) X z 

r- ( 1 ~ X i 
-Li 2 

-Li 2 



Li 2 



1 - A, 



1 - Xi - X z p~ 
(1 - A,) 2 



Li 2 



1- Xi - X z p+ J 
' \-Xi-Xz 



Li 5 



' \-\j-\z 
\-Xi - XzP-i 



I - X z - X z p- 



^((i-X^l-h-Xz) 



{1-Xif + XiXz) V (l-A,) 2 + A,Az 
i yi dxdyy 2 [l - yx(l - x)) 1 [5 2(1 - A 



o Jo 
Xi 



BtiXi) 
2A 



^-LlnAi-- ^lnAi + 4 

Az 1 — Aj \ 

2(1 -A, -A Z )(1-A,) + A,A; 
X z 



X, 
'1 - A, 



A; 



A; 



1 r\ tan 



-i 



f I dxdy y In B^Xi) = - 
'0 Jo 2, 



Xi (1 — Aj) + XiX z T(x v 

-—In Aj H I(Xi), 

Az A^ 



3 1-Aj / 2(1 -Xi 



(A.2) 



(A.3) 



r\ tan 1 



(A.4) 



Ki(Xi, Xj 



1 /•! 



c?X(iy y 1 
o 5 2 (Ai,Aj) X z 
1 - A,- \ 



+Li 2 



1 — Aj + X Z ^ 



Li 2 
Li 2 



1-Aj 



^1-Aj + A z ^ + 
1-Aj + A z 
1 - A, + X z ^ 



( l-Xj+Xz 
Ll2 U-A. + A^ 
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-Lb 



(l_A.)(l_ Aj .) + A Zi 



(l_A.)(l_ Aj .) + Az, 



^ 2 (A,,A,) = I' f 
Jo Jo 



dxdy y 2 
o Jo B 2 (Xi,Xj) 
1 A,; — A, 

2 



1 

Az 

A, 



- 1 + 



2A Z y 
+(l-A j )K 1 (A i ,A J ) 



ln(^| - 



A, / Xz 



1-A, 



tan 



In A ?: 



Aj + A, — A^ 



+ (Aj <-> A 



K(A,,A,) = /V 



1 Z" 1 dxdyy 3 x(l — x) 1 

~A^ 



E>2(Xi, Xj) 



3)i 



1 2 — A, — Ao 
- H 

2 Az 



1 , M (2-A,-A J - + Az)(A 3 --AQ /V 
A^ ln(AlAj) + 2A| ln \X 3/ 

2-A,-A jv ^ tan - 1 



A; 



A; 



Ai + A, — A; 



^2(Aj,Aj) — 



dxdy y in B 2 {Xi, Xj) = -- - ■ 



(A.5) 



(A.6) 



(A.7) 



io Jo 

Aj - A,- 



_^ + , (i + __ )ln(A(Aj) 



+ 



A,; 



4A| 



J (2 + 2A z -A J -A J )lnl ^] 4 



2 - Aj - Xj + A ; 



2A; 



<w _-, 
x-t — tan 

X7, 



w 



. Ai + Xj — X2 



+ 



. '3/ "»Z 

(1- Aj)(l- Aj) + A z 
Az 



^(A^A,), (A.8) 



where p± — (1 ± 277) /2 with 77 = y 4A Z 1 — 1, £± = (Az — Aj + Aj ± i^/w)/2Xz with u> 
4AiAj — (Az — Aj — Aj) 2 , and 



B^Xi) = (l-y)X l + y[l-X z yx(l-x)], 
B 2 (Xi,Xj) = l-y + y[xXi + (l-x)Xj-X z yx(l-x)]. 



(A.9) 
(A.10) 



Note that w > for |\Ai — yXj\ < vAz < vA7+ y Aj. If vA7+ y Aj < yAzi then one has 
to analytically continue the function 



w tan 



A, + Aj — A^ 



2\/wtan 1 



X z ~ (VA7 - JXj) 



\l (Vxl + ^XjY-Xz 



a/— w ln 



Az 



Az 



, _ Aj) 2 
1 A z 



(v^+AAi) 2 



— m^/— 



(A.ll) 



Az 
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The dilogarithmic function Li2(x) (with x being real) should also be continued analytically 
as follows: 

U 2 {x±ie) = - f dt—^— 4 ± i9(x- l)7rlnar. (A.12) 
Jo t 



For VAz < | vAl — \/A?Ij we nave checked that Eq. (A. 5) and the l.h.s of Eq. (A. 11) do 
not contain any imaginary part. This implies that JF|' is pure dispersive in this specific 
kinematic range. 

As we are interested in heavy neutrinos with masses larger than Mz, the absorptive 
part of T l z will solely originate from the Fig. 2(i) in which only intermediate light neutrinos 
can come kinematically on-mass shell. Neglecting light neutrino masses in the calculation, 

we get 



Abs(J?) 



i=l 



3 
2 



1 



1 + — ) ln(l + A z ) 



(A.13) 



B Loop functions of flavour-violating decays 
of charged leptons 



In Section 4 the amplitudes of the flavour-violating decays of I, I — ► Vl\li and I — > l ,r y, 
are expressed in terms of all possible form factors that are derived by an explicit calculation 
of the Feynman graphs shown in Figs. 1-3. The photonic form factors, F 1 and G 1 in 
Eq. (4.1), vanish in the limit of zero external momenta and lepton masses due to the 
electromagnetic gauge invariance. One has consistently to expand the corresponding loop 
integrals up to the next order of q 2 [57| in order to obtain a nonvanishing result. After a 
straightforward computation, we find that 



F 7 (x) 



Gr 7 (x) 



F z (x) 
G z (x,y) 
Hz(x,y) 



7x 3 



x 



12x 



x 



10x 3 + 12x 2 



x 



12(1 



2x 3 + 5x 2 



x 



6(1 

3x 3 



x 



lnx, 



4(1 

5x 



x 



2(1 



x 



2(x - y) 
Vxy 



x) 3 2(1 

5x 2 

2(1 -x) 2 

x 2 {l - y) y 2 
mx — 



Inx, 



lnx, 



Box 



[x,y) 



4(x - y) 
1 



1 — x 

x 2 — Ax y 1 

1m- 

1 — x 



1 - 
Ay 



x-y 



xy . 
1 + f) 



1 — x 



+ 



1-2/ 
x 2 hax 
(1 - x)' 



-y 

In y 



In y 



_ _ y 2 Iny 
y (i-y) 2 



(B.l) 
(B.2) 
(B.3) 
(B.4) 
(B.5) 
(B.6) 
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( 1 x lnx 

-2xy I h 
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1 — x (1 — x) 2 I — y (1 — y) 2 
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(4 + xy) 
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v l-x (1-x) 2 1-y (1-y) 2 , 
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1-x ' (1-x) 2 1-y (1-y) 2 / J' ^'^ 

Although F 7 , G 7 , F^, and -Fb 02 , are already known in the literature p7H55[p5[[ , the form 
factors Gz, H z and Gbox are newly obtained by Eqs. (B.4), (B.5) and (B.7), respectively. 

For completeness, we list below expressions of the form factors computed at some 
special values of the arguments 
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(B.8) 
(B.9) 
(B.10) 
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(B.ll) 
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0, # z (0,0) = 0, # z (l,0) = 0, If z (l,l) = ~; (B.12) 
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(B.13) 
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G Box (l,x) 



—\/x 



x 3 - 2x 2 + 7x - 6 x 2 - 4x , 

— h — — mx 



2(1 -xf 



(1-xf 



G Box (0,x) = 0, G Box {0,0)=0, G Box (l,0) = 0, G Boa: (l, 1) = -. (B.14) 

Since all the form factors given in Eqs. (B.1)-(B.7) are multiplied by certain combi- 
nations of B and C matrices in the decay amplitudes (4.1), (4.6), (4.7) and (4.8), it will 
be helpful to define the following composite form factors: 



F" = Y B ii B i'i F i(^i) = J2Bi N .Bi /Ni F y (X Ni ), 
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(B.18) 



where we have made use of the identities of Eq. (2.9) in the final step of the Eqs. (B.15)- 
(B.18) and re-expressed all the composite form factors as a sum over the heavy neutrino 
states. This simplification enables us to study the behaviour of these form factors in the 
heavy neutrino limit. 

For the purpose of illustration, we will discuss the results of this asymptotic limit in a 
model with two heavy Majorana neutrinos. Employing Eqs. (5.1) and (5.2) for the mixing 
matrices B and C, we find that for A^ = m^/M^ 3> 1 and p = m 2 N2 /m 2 Ni 1, 
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t Box ~ ( S L S L °hh + S L S L Ol'h) 




(B.22) 



In the limit p — > 1 and for 3> 1, Eqs. (B.21) and (B.22) take the form 




3 1 v n ° 



(B.23) 



i=l 



Box 



(B.24) 



From Eqs. (B.19)-(B.24), it is obvious that all the composite form factors, F3 , G l J: , F§ , 
and F^^ 2 , violate the decoupling theorem pg| . Note that terms of C((s^) 2 ) in Ff' depend 
logarithmically on the heavy neutrino mass, m^, while terms proportional to (s^) 4 in 
Eqs. (B.23) and (B.24) show a strong quadratic, m% i /M 2 v , dependence and should not be 
neglected in the calculation when m^ 1 > 200 GeV (see, e.g., Fig. 4). 

C Three-body phase-space integrals 

As we have seen from Eq. (4.6), the 7-mediated amplitude of the decay l(p) — > 
I' {p')h{Pi)h{P2) contains a non-local interaction which leads to a collinear singularity in 
the limit q 2 = (pi + P2) 2 — ► 0. This divergency can only be avoided if one assumes 
that the leptons, l\ and 1%, coupled to the virtual photon are not strictly massless, i.e. 
mi x = mi 2 = e ^ 0. Thus, after performing phase-space integration, we neglect all those 
terms that vanish as e — > 0. On the other hand, the mass of V can safely be set to zero. 
Then, the phase-space boundaries can be given by 



where si = (p' + P2) 2 , $2 = (pi + P2) 2 ■, m is the mass of the decaying lepton /, and s 1 is 
the upper (lower) limit of the Mandelstam variable s\. 

The divergent phase-space integrals relevant for the decay I — > VI1I2 (with I' 7^ I2) are 
given by the following expressions: 





(C.l) 



Pi 




(C.2) 
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P 2 = J J ds2dsi s J. = m ^^i n ^-- 7 -) + (e), (C.3) 

P 3 = f f ds2dsi i = m \]-\n^-^) + 0(e). (C.4) 

Note that a different result would have been obtained in Eqs. (C.2)-(C.4) if we had origi- 
nally expanded the square root existing in sf in terms of e and then performed the phase- 
space integration. Apparently, this technical problem seems to have caused some confusion 
in the literature, as far as the correct analytic expression of the non-local interaction in 
Eqs. (4.9) and (4.10) is concerned. In the three-body leptonic decays of / where I' = l\ = I2, 
one may have to take into account an additional divergent phase space integral when in- 
terfering the two possible, si-channel and s 2 -channel, 7-exchange amplitudes, i.e. 

r r \ jr^ 1 Tfi^ 1 

P 4 = ds 2 d Sl = ---ln 2 2 + -ln 2 — + 0(e). (C.5) 

J J S2S1 ±2 Z £ 

The integral P 4 in Eq. (C.5), however, is always multiplied by the small mass e of the final 
leptons and therefore goes to zero as e — > 0. As a result, the only type of divergency for 
e — > that appears in Eqs. (4.9) and (4.10) is the logarithmic one, \n(m 2 /e 2 ). 
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Figure Captions 

Feynman graphs responsible for generating the effective vertex jll' (I 7^ V). 

Feynman graphs responsible for generating the effective vertex ZIV (I 7^ V) . 

Feynman diagrams relevant for the leptonic decays I — > I'l^- 

B(t~ — > e~e~e + ) (solid line) and B(t~ — > e~ fi~fi + ) (dashed line) as a func- 
tion of the heavy neutrino mass mjv(= = itin 2 ) assuming {s v £) 2 = 0.07, 
( s l) 2 = 0.015 and (s^ 1 ) 2 ~ 0. Numerical results obtained when seem- 
ingly suppressed terms of C((s^) 4 ) are neglected, are also presented for 



e pi /i + ) (dash-dotted line). 
e~ fi~fi + ) (dashed line) as a func- 



B{t — > e e e + ) (dotted line) and 5(r 

Z?(t~ — > e~e~e + ) (solid line) and .B(t _ 
tion of the heavy neutrino mass mjv(= = m^) using (s^) 2 = 0.035, 
( s l) 2 = 0.010 and (s^ M ) 2 ~ 0. We also display numerical results obtained 
by neglecting seemingly suppressed terms of (9((s^) 4 ) in the calculation of 
B(t~ 

B(t~ — > e e e 

(solid line) and 500 GeV (dashed line). We have assumed (s£ T ) 2 = 0.07, 
( s l ) 2 = 0.015 and (s^ 1 ) 2 — 0. The corresponding numerical results of 
B{t~ — > e~/i~ii + ) are shown for = 200 GeV (dotted line) and 500 GeV 
(dash-dotted line). 

Numerical estimates of B(Z — > e~r + ) + S(Z — > e + r~) as a function of the 
heavy neutrino mass m^{— = m^ 2 ) for three representative values of 



e e e + ) (dotted line) and B{r — > e /i yU + ) (dash-dotted line). 

+ ^ as a function of the ratio m^ 2 /m^ 1 for — 200 GeV 



the mixing parameters (sx, T ) 2 and (s£ 
(s£) 2 = 0.015 (solid line), (ii 
line), and (iii) (s u £f = 0.020 and (s 



L ) 



0): 



0.070 and 



^) 2 



0.035 and (s£) 2 = 0.010 (dashed 
0.010 (dotted line). 

Numerical estimates of B[Z — > e~r + ) + B(Z — > e + r~) versus m^/m^ 
for selected values of = 200 GeV (solid line), 400 GeV (dashed line), 
600 GeV (dotted line), and 1 TeV (dash-dotted line). We have used (s v £ f = 
0.07 and (s^) 2 = 0.015. 
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